A1 = readmatrix ('NiI2_CPGE_Spectra_-24V_30K_600-650_2.csv')
A2 = readmatrix ('NiI2_CPGE_Spectra_-24V_30K_650-700.csv')
A3 = readmatrix ('NiI2_CPGE_Spectra_-24V_30K_700-750.csv')
%B = readmatrix('supercontinuum_power_calibration.csv')

x1 = A1(:,1)
y1 = A1(:,3)*1e12*4.3
y2 = A1(:,5)*1e12*4.3
y3 = A1(:,6)*1e12*4.3

x2 = A2(:,1)
m1 = A2(:,3)*1e12
m2 = A2(:,5)*1e12
m3 = A2(:,6)*1e12

x3 = A3(:,1)
n1 = A3(:,3)*1e12*0.7
n2 = A3(:,5)*1e12*0.7
n3 = A3(:,6)*1e12*0.7

hold on
plot(x1,y1,LineWidth=1,Color=[0.8 0 0])
plot(x1,y2,LineWidth=1,Color=[0 0 0.8])
plot(x2,m1,LineWidth=1,Color=[0.8 0 0])
plot(x2,m2,LineWidth=1,Color=[0 0 0.8])
plot(x3,n1,LineWidth=1,Color=[0.8 0 0])
plot(x3,n2,LineWidth=1,Color=[0 0 0.8])
%plot(x,z1,LineWidth=1,Color=[0.8 0 0])
%plot(x,z2,LineWidth=1,Color=[0 0 0.8])
%plot(x,n,LineWidth=1,Color=[0 0 0.8])
%plot(x,r,LineWidth=1,Color=[0 0 0])
xlim([600,750])
xticks(600:50:750)
%ylim([0,0.11])
%yticks([0:0.05:0.1])
set(gca,'XMinorTick','on','YMinorTick','on')
set(gca,'box','on','Xcolor',[0 0 0],'Ycolor',[0 0 0])
hold off
box on
set(gcf,'Position',[0, 0, 200, 200])